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Iy  =  Ox  *  V,  -  0„  *  V, 

Ij  =  O*  *  Vy  -  Oy  *  Vt. 

The  inclination  angle,  INCL,  of  the  SPAS  is  then 


INCL  -  cos'1 
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The  SPAS  heading  angle  measured  from  true  north,  HEAD,  is  computed  from: 


HEAD 


tan'1 


_ cos ( INCL) _ 

[sin2  {INCL)  -  sin2  (GEOCLAT)  j  1/2 


* 


If  the  SPAS  is  traveling  from  north  to  south,  determined  by  examining  the  latitudes,  the 
supplement  of  the  heading  angle  is  used.  The  azimuth  from  true  north  and  elevation  measured 
from  the  local  horizon  of  the  SPAS  line-of-sight  are  determined  as  follows: 

AZIMUTH  =  YAW  +  HEAD 

ELEVATION  =  PITCH. 

The  Julian  date,  JD,  is  computed  from: 

JD  =  INT  (365.25  y)  +  INT  (30.6001  (m+1))  +  DD.dd  +  1720994.5+  B 

where  INT  is  the  integer  part  of  the  value  in  parentheses  and  DD.dd  is  the  decimal  day,  y  is  the 
year  and  m  is  the  month  at  the  time  of  the  calculation. 

B  =  2  -  A  +  1 4  I  /  ancf 

4 


The  time  T,  measured  in  Julian  Centuries  of  36525  ephemeris  days  from  the  epoch  1900 
January  0.5  ET,  is  given  by 


T  -  <777-  2415020.0 
36525 


Since  T  is  expressed  in  centuries,  an  error  in  the  5th  decimal  place  will  cause  an  error  of  0.4 
days  in  the  time.  Therefore,  it  is  essential  that  this  calculation  be  done  in  double  precision. 

v  J  * 

The  geometries  mean  longitude  of  the  sun,  L,  referred  to  as  the  mean  equinox  of  the  date, 
is  given  by:  « 


L  ~  p 79.69668  +  36000.76892  T  ♦  0.0003025  T2 . 


•  The  mjjan  anomaly  of  the  sun,  M,  is  expressed  as: 


M  =  358.47583  +  35999.04975  T  -  0.000150  T2  -  0.0000033  T3. 
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The  equation  for  the  center  of  the  sun,  C,  is  determined  from: 

C  =  (1.919460  -  0.004789  T  -  0.000014  T2)  sin(M)  +  (0.0200943  - 
0.000100  T)  sin(2M)  +  0.000293  sin(3M). 

The  true  longitude  of  the  sun,  0,  equals 

0  =  L  +  C. 

This  longtitude  must  be  modified  to  obtain  the  apparent  longitude  of  the  sun,  referred  to 
as  the  true  equinox  of  the  date.  0  is  corrected  for  mutation  and  aberration  as  follows: 

Q  =  259.18  -  1934.142  T 
0^  *  0  -  0.000569  -  0.00479  sin  (0)  . 


The  obliquity  of  the  ecliptic,  e,  is  given  by 


e  =  23.452294  -  0.0130125  T  -  0.00000164  T 2  + 
(9.2100  ♦  0.00091  T)  COS  (Ql/3600.0. 


The  right  ascension  at  the  sun,  y,  is  given  by  the  formula 

|  (hours)  . 

p  p'  ) 

The  declination  of  the  sun,  6,  is  given  by 
sin  (6)  =  sin  (0^)  sin  (e). 

The  coordinates  of  the  sun  relative  to  the  position  of  the  satellite  must  be  determined. 

The  sidereal  time  at  Greenwich  at  0  hr.  UT  expressed  in  decimal  hours  in  the  interval  0-24 
hours  is  given  by  the  formula: 

TIME*  =  0.276919398  +  100.0021359  T. 

The  sidereal  time  at  Greenwich  at  the  time  of  the  calculation,  Ts  is  therefore: 

Ts  cal  =  1.00273908  *  TIME  rA1.  +  TIMEs 


Y  = 


180 

15it 


tan 


cos  e  sm 


cos  (0 
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where  TTME,^  is  in  decimal  hours. 
The  local  hour  angle  in  degrees  is 


LHA  =  Tscal  -  SPAS  LON  -  y. 

The  angle  A  in  Fig.  I,  which  is  the  compliment  of  the  solar  zenith  angle,  is  given  by: 
A  =  tin'1  [sin  (GEOCLAT)  sin  (6)  +  cos  (GEOCLAT)  cos  (6)  cos  (LHA)]. 


Figure  1.  SUN -SPAS -EARTH  Geometry. 


0 


* 


The  earth  center  angle  between  the  SPAS  and  the  sun,  the  solar  zenith  angle  of  the  SPAS, 
is 

SUNj^  =  90*  -  A 


The  azimuth  to  the  specular  point,  i.e.,  azimuth  of  sun,  is  defined  as: 


17  ■  fan'l  |  Sin  {LHA)  I  +  ■) 

SPEC  ‘COS  (LHA)  sin  (GEOCLAT)  -  tan  (b)  cos  ( GEOCLAT ) 


Five  specular  point  altitudes  are  assumed,  specifically  40,  30,  20,  10  and  0  km.  An 
iterative  procedure  is  necessary  to  calculate  the  earth  center  angle  between  the  SPAS  and  the 
specular  point. 
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With  B  equal  to  SUN^.  HsrerL1AR  equal  to  one  of  the  5  specular  point  altitudes  measured 
from  the  center  of  the  earth,  and  B'  initially  equal  to  B  the  following  parameters  are  computed: 

Xx  -  2  B’  -  B 

FN  =  Hspec-ula-±  sin  <B')  -  sinUj) 

^ spas 

FN '  =  -^^ESSilSl  cos  (S')  -  2 COS  (Xx)  . 

^  spas 


A  check  is  made  to  see  if  (B-B')  is  less  than  0.0001.  If  it  is  not  X„  FN  and  FN'  are 
recomputed  using  the  new  B'.  This  process  is  repeated  until  (B^,  -  B'n)  becomes  less  than  0.0001. 
The  earth  center  angle  between  the  satellite  and  the  specular  point  is  then  given  by: 

ECAsAX.sreiC  =  B  -  B  . 

The  formalae  used  in  this  iteration  can  result  in  a  B'  that  is  greater  than  360°.  Therefore, 
B'  must  be  reduced  by  multiples  of  360°  until  an  angle  less  than  360°  is  reached. 

Figure  2  illustrates  the  geometry  described  above. 

The  elevation  to  the  specular  point  is  given  by: 

ELEVspbc  =  -[90*  +  SUN^  -  2  B']. 

If  B'  is  greater  than  90°,  there  is  no  specular  point  for  the  given  input  parameters. 

The  distance  D  in  kilometers  from  the  SPAS  to  the  specular  point  is  given  by: 

D2  ~  HsPAs  +  ^SPECULAR  ~  ^  ^SPAS  ^ SPECULAR  CC>S  (ECA gAT-SPBC ^  ’ 
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-  B'-X, 

-8-B' 


Figure  2.  Specular  Scattering  Geometry 
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The  latitude,  LATsrec,  of  the  specular  point  is  defined  as: 

LATSPSC  =  sin'1  [cos  ( ECAsat_ spgc)  sin  (GEOCLAT)  * 
sin  (ECAsat_spbc)  cos  (GEOCLAT)  cos  (AZspbc)  ]  . 


The  longitude  of  the  specular  point  is  determined  in  the  following  manner.  First  the 
parameters  A,  and  A,  are  computed  from: 


Ax 


sin"1 


sin  (  ECAsat_spbc)  siti  (AZSPSC) 
COS  ( latspbc ) 


cos  ( ECAsat_spbc)  -  sin  (GEOCLAT)  sin  { LATspec ) 
COS  < GEOCLAT)  COS  ( LATspbc)  . 


Aj  is  the  cos  of  A„  however  only  the  sign  of  A2  is  required  to  determine  X,  the  longitudinal 
difference  LONGSPAS  -  longsrec. 

If  A,  0°  and  As  <  0,  then  X  =  7t  -  A, 

If  At  <  (f  and  Aj  <  0,  then  X  =  -n  -A,. 

If  A^  >  0,  then  X  =  A,. 

The  longitude  of  the  specular  point  is  then 
LONGsrec  =  SPASLON  -  X. 

Since  LONGsrec  must  be  between  0  and  36 0°,  the  same  modification  must  be  made  to  it 
as  was  done  for  B'  to  insure  that  it  stays  within  this  range. 

The  pointing  angle  to  the  specular  point  is  defined  as: 

P  =  cos1  [cos  (ELEVspeo)  cos  (ELEVATION)  cos  ( AZIMUTH- AZs,^)  + 
sin  (ELEVsrec)  sin  (ELEVATION)] 
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The  results  of  a  typical  calculation  are  shown  in  Table  1.  The  first  column  contains  the 
day  of  the  year  (0  -  365),  the  hours  of  the  day,  the  minutes  of  the  hour  and  the  seconds  of  the 
minute.  The  time  is  given  in  GMT.  The  next  column  contains  the  mission  elapsed  time  (sec) 
The  third  column  which  contains  the  elapsed  time  (sec)  from  the  start  of  an  event  was  included 
for  plotting  purposes.  The  last  5  columns  contain  the  pointing  angles  from  the  SPAS  position  at 
the  time  shown  in  column  1  to  a  specular  point  seen  at  40,  30,  20,  10  and  0  km.  The  number 
9999.9  in  one  of  these  columns  indicates  that  the  SPAS  cannot  see  a  specular  point  at  that  altitude. 
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TABLE  1 


EXAMPLE  OF  OUTPUT 


MISSION 

HEIGHT  OF 

THE  SPECULAR  POINT  1 

ELAPSED 

ELAPSED 

DAY  HH : MM: SS 

TIME  (SEC) 

TIME 

40  KM 

30  KM 

20  KM 

10  KM 

0  KM 

121/23:42:11,00 

302937.00 

r  ' 

66 

14.95 

15.01 

15.08 

15.14 

15.20 

121/23:42:12.00 

302938.00 

L  . 

14.90 

14.97 

15.03 

15.09 

15.15 

121/23:42:13.00 

302939.00 

68 

14.85 

14.92 

14.98 

15.04 

15.11 

121/23:42:14.00 

302940.00 

69 

14.81 

14.87 

14 . 93 

14 . 99 

15.06 

121/23:42:15.00 

302941.00 

70 

14.76 

14.82 

14.88 

14  .  95 

15.01 

121/23:42:16.00 

302942.00 

71 

14.71 

14.77 

14.84 

14 . 90 

14 . 96 

121/23:42:17.00 

302943.00 

72 

14.66 

14.73 

14.79 

14 . 85 

14.92 

121/23:42:18.00 

302944.00 

73 

1 

14.68 

14.74 

14.80 

14 . 87 

121/23:42:19.00 

302945.00 

74 

EBB 

14.63 

14 . 69 

14.76 

14.82 

121/23:42:21.00 

302947.00 

wm 

14.48 

14.54 

124.60 

14.67 

14.73 

121/23:42:22.00 

302948.00 

77 

14.43 

14.49 

14.56 

14.68 

121/23:42:23.00 

302949.00 

78 

14.38 

14.44 

14.51 

1 

14.64 

121/23:42:24.00 

302950.00 

mm 

14.33 

14.40 

14.46 

14 . 52 

14.59 

121/23:42:25.00 

302951.00 

80 

14.28 

14.35 

14.41 

14.48 

14 . 54 

121/23:42:26.00 

302952.00 

81 

14.24 

14.30 

14.36 

14.43 

14.49 

121/23:42:27.00 

302953.00 

82 

14.19 

14.25 

14.32 

14.38 

14.45 

121/23:42:28.00 

302954 . 00 
_ _ 

83 

14.15 

14.21 

I  '  '  1 

14.34 

14.41 

121/23:42:29.00 

302955.00 

84 

14.10 

14.17 

mm 

14.30 

14.36 

121/23:42:30.00 

302956.00 

85 

14.06 

14  . 12 

14.19 

14.25 

14.32 

121/23:42:31.00 

_ 

86 

14.01 

14.07 

14.14 

- .  j 

BEB 

14.27 

121/23:42:32.00 

302958.00 

87 

13.96 

14.03 

14.09 

14.16 

14.22 

121/23:42:33.00 

302959.00 

88 

13.92 

13.98 

14.05 

14.11 

14 . 18 

121/23:42:34.00 

302960.00 

89 

13.87 

13.94 

14.00 

14.07 

14.13 
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